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Abstract 

Mixed-spin Ising model on a decorated Bethe lattice is rigorously solved by 
combining the decoration-iteration transformation with the method of exact re- 
cursion relations. Exact results for critical lines, compensation temperatures, 
total and sublattice magnetizations are obtained from a precise mapping rela- 
tionship with the corresponding spin-1/2 Ising model on a simple (undecorated) 
Bethe lattice. The effect of next-nearest-neighbour interaction and single-ion 
anisotropy on magnetic properties of the ferrimagnetic model is investigated in 
particular. It is shown that the total magnetization may exhibit multicompen- 
sation phenomenon and the critical temperature vs. the single-ion anisotropy 
dependence basically changes with the coordination number of the underlying 
Bethe lattice. The possibility of observing reentrant phase transitions is related 
to a high enough coordination number of the underlying Bethe lattice. 
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1. Introduction 

Fcrrimagnetism has been extensively studied in the past both experimen- 
tally as well as theoretically. The magnetic structure of insulating ferrimagnets 
is built up from several (at least two) sublattices, which are occupied by different 
carriers of unequal magnetic moments interacting with each other through the 
antifcrromagnetic interaction. Mixed-spin Ising models probably represent the 
simplest theoretical models, which may bring an insight into the most important 
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features accompanying ferrimagnetic ordering such as for instance a compensa- 
tion phenomenon. A compensation point denotes the temperature T compi at 
which the resultant magnetization vanishes due to a complete cancellation of 
sublattice magnetizations below the critical temperature T c , i.e. T comp < T c 
[1]. An existence of compensation temperatures has obvious technological sig- 
nificance in the field of thermomagnetic recording, because only a small driving 
field is needed at a compensation point in order to achieve a magnetic pole 
reversal @, Q. Besides, the mixed-spin Ising models belong to the most inter- 
esting extensions of the standard spin- 1/2 Ising model, which may display more 
diverse critical behaviour compared with their single-spin counterparts. 

Decorated mixed-spin Ising models, which were originally introduced by 
Syozi and co-workers (see Ref. [4J and references cited therein), have been inten- 
sively studied many years ago as theoretical models of ferrimagnetism. Thermal 
variations of the magnetization in the decorated Ising ferrimagnets have been 
investigated in detail and these results shed light on some characteristic fea- 
tures of real ferrimagnetic materials Q . The effect of single- ion anisotropy 
on magnetic properties of the mixed spin-1/2 and spin-S 1 (S > 1/2) Ising ferri- 
magnet on decorated planar lattices has been examined by Kaneyoshi with the 
help of differential operator technique and the effective- field theory 0, H| . By 
making use of these theoretical tools, Kaneyoshi has predicted many interesting 
features of the decorated Ising ferrimagnets such as an appearance of two com- 
pensation points and/or multiple reentrant phase transitions. It is worthwhile to 
remark that Kaneyoshi's predictions were subsequently confirmed by Jascur Q 
and Dakhama [13] when solving the ferrimagnetic version of the decorated Ising 
model rigorously. The significance of this model also lies in the fact that tem- 
perature dependences of the total magnetization show numerous characteristic 
features not predicted in the standard Neel theory of ferrimagnetism [l[. 

Recently, the magnetic properties of various decorated mixed-spin Ising 
models have been explored in detail by variety of mathematical techniques 
with some exact K, 11, 12, 13, 14, 15,(3] as well as approximate results 

26 1 . The present work deals with magnetic 
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properties of the mixed-spin Ising ferrimagnet on the decorated Bethe lattice 
accounting for the next-nearest-neighbor interaction and the uniaxial single-ion 
anisotropy. The exact solution will be obtained by combining the decoration- 
iteration transformation with the another rigorous method based on recursion 
relations. First, the decoration-iteration transformation will be applied to es- 
tablish the precise mapping correspondence between the mixed-spin Ising model 
on the decorated Bethe lattice and its equivalent spin-1/2 Ising model on the 
simple (undecorated) Bethe lattice. It is well known that this latter model 
can be subsequently treated by making use of recursion relations [27j ]. which 
will help us to complete our exact calculation for the original mixed-spin Ising 
model on the decorated Bethe lattice. Within the framework of these two ex- 
act methods, we will focus on the effect of the single-ion anisotropy and the 
next-nearest-neighbour interaction on the phase diagrams, critical behaviour 
and compensation phenomenon of the model under investigation. 

The outline of the present paper is as follows. In Section 2, the detailed 
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description of the model system will be presented and the basic steps of both 
exact methods will be briefly clarified. The most interesting results are presented 
and discussed in detail in Section 3. In particular, our attention is focused on the 
finite-temperature phase diagrams, reentrant phase transitions and temperature 
variations of both sublattice as well as total magnetizations. Finally, some 
conclusions are mentioned in the last Section IV. 



2. Model and method 

Let us begin by considering the mixed spin-1/2 and spin-S 1 Ising model on the 
decorated Bethe lattice with a quite general coordination number q as schemat- 
ically illustrated in Fig. Q] on the particular example of the Bethe lattice with 
the coordination number q = 3. The magnetic structure of the investigated 
model system constitutes the Ising spins a — 1/2 placed on lattice sites of a 
deep interior of infinite Cayley tree (Bethe lattice) and the Ising spins of a quite 
general magnitude S placed on each bond of the original Bethe lattice. The 
total Hamiltonian of the mixed-spin Ising model on the decorated Bethe lattice 
then reads 

nn nnn Nq/2 

u = - J E s «>i - f E a w - D E s l « 

i,j k,j i=l 

where aj = =Fl/2 and Si — —S, —S + 1,...,S, the relevant subscript specifies 
the lattice position and the exchange constants J and J' represent the nearest- 
neighbour interaction between the spin-1/2 and spin-S' atoms and, respectively, 
the next-nearest-neighbour interaction between the spin-1/2 atoms. Further, 
the parameter D stands for the uniaxial single-ion anisotropy acting on the 
spin-5 atoms only. As one can see from Fig. [TJ the magnetic structure of the 
investigated model system constitutes two interpenetrating sublattices. The 
former sublattice A is formed by the sites of original (undecorated) Bethe lattice 
that is occupied by the atoms with the fixed spin a = 1/2, while the latter 
sublattice B is occupied by the decorating atoms with an arbitrary spin value 



S. The calculation on the Bethe lattice is done recursively [27 1. 

Now, we will turn our attention to main points of the transformation method, 
which provides a rigorous mapping equivalence between the investigated spin 
model and its corresponding spin-1/2 Ising model on undecorated Bethe lattice. 
First, it is quite useful to rewrite the total Hamiltonian (JXJ) as a sum of the bond 
Hamiltonians, i.e. H — J2^=( 2 ^k, where N denotes the total number of the 
spin-1/2 atoms from the sublattice A, q is the coordination number of the un- 
decorated Bethe lattice and the summation runs over all bonds of the decorated 
Bethe lattice. Thus, each bond Hamiltonian Hk involves all the interaction 
terms including the decorating Ising spin Sk and the next-nearest-neighbour 
interaction between its two nearest-neighbouring Ising spins Uki and a k 2 (for 
better illustration see spin cluster enclosed by ellipse in Fig. [TJ) 

Hk = —JSk(cTkl + a k2) - DS k — J'<TklVk2- (2) 
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Figure 1: The mixed spin-1/2 and spin-S 1 Ising model on the decorated Bethe lattice (figure on 
the left) and its relation to the corresponding spin-1/2 Ising model on the simple (undecorated) 
Bethe lattice with the coordination number q = 3 (figure on the right). The blue (red) circles 
denote lattice positions of the spin-1/2 (spin-1) Ising spins and the ellipse demarcates all the 
interaction terms belonging to the fcth bond Hamiltonian. 



A crucial step in our procedure represents the evaluation of the partition func- 
tion. For further convenience, the partition function of the model under consid- 
eration can be written in the following useful form 

Nq/2 

2 = EIlE ex pW' ( 3 ) 

{cTi} k=l S k 

where (3 — l/(fcsT), ks is the Boltzmann's constant and T is the absolute 
temperature. The symbol J2s k denotes a summation over 25 + 1 spin states 
of the kth decorating Ising spin from the sublattice B, while the symbol } 
marks a summation over all available spin configurations of the complete set 
of the spin-1/2 Ising atoms from the sublattice A. It is noteworthy that the 
summation over spin states of the decorating spin-5 atoms can be performed 
independently of each other because of the lack of any direct interaction among 
them. After performing the summation over the spin states of the decorating 
Ising spin Sk, one gets the effective Boltzmann's weight that implies a possibility 



of applying the generalized decoration-iteration transformation [281 . 1291 . |3 ( 



^exp(-/3?4) = exp(/3JV fe i<7 fc2 ) ^ exp((3Dn 2 ) cosh[(3 J n(cr kl + a k2 )) 

S k n=-S 

= Aexp(PRa kl a k2 ). (4) 

The physical meaning of the mapping transformation Q lies in removing all 
the interaction parameters associated with the decorating spin S k and replac- 
ing them by a new unique effective interaction R between its two nearest- 
neighbouring Ising spins o k \ and a k 2- It is worthwhile to mention that the 
both mapping parameters A and R are 'self-consistently' given directly by the 
transformation formula (J4|, which must hold for all four possible spin combina- 
tions of two Ising spins a k \ and o k i that provide just two independent equations 
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from the mapping transformation (U). In the consequence of that, two yet un- 
known mapping parameters A and R can be unambiguously determined by the 
following formulas 

A=(V!V 2 )i, pR=l3J' + 21n(^f), (5) 



which contain two newly defined functions 

S s 
Vi = exp((3Dn 2 )cosh((3Jn), V 2 = ^ exp(/3Dn 2 ). (6) 

n—~S n——S 

After straightforward substitution of the mapping transformation Q into 
the formula ([3]) , one easily obtains an exact relation between the partition func- 
tion of the mixed-spin Ising model on the decorated Bethe lattice and respec- 
tively, the partition function of the spin-1/2 Ising model on the simple (un- 
decorated) Bethe lattice with the temperature-dependent effective interaction 
R 

Z{fi, J, J', A S, q) = A N «/ 2 Z Bethe (0R). (7) 

The mapping relation ([7]) between both partition functions represents a central 
result of our calculation, because this relationship can in turn be employed for 
a rigorous calculation of some important quantities such as sublattice magneti- 
zations and free energy. It is worthwhile to recall that this mapping correspon- 
dence is valid for arbitrary values of the decorating spin S and also independently 
of the lattice coordination number q. 



Exact mapping theorems developed by Barry et al. [311 . 132t 1331 . [34j in com- 
bination with Eq. ([7]) yield a straightforward relation between the spontaneous 
magnetization of the sublattice A of the mixed-spin Ising model on the dec- 
orated Bethe lattice and respectively, the spontaneous magnetization of the 
corresponding spin-1/2 Ising model on the undecorated Bethe lattice 

m A = i ((erfei) + (<r k2 )) = i ((er fe i) + (o-fc 2 )) B cthe = m Bc thc(/3i?)- (8) 

Here, the symbols (• • •) and (• • •) Bethe mean the standard canonical ensemble 
average performed in the mixed-spin Ising model on the decorated Bethe lattice 
defined through the Hamiltonian ((TJ) and respectively, its corresponding spin- 
1/2 Ising model on the undecorated Bethe lattice. According to Eq. ([S]), the 
sublattice magnetization tua directly equals to the spontaneous magnetization 
?TiBcthe of the corresponding spin-1/2 Ising model on the simple Bethe lattice 
provided that Eqs. JS])-® lead to the effective ferromagnetic interaction R > 0. 
In the case of the effective antiferromagnetic interaction R < 0, one alternatively 
gains a quite similar exact result for the staggered magnetization m s A of the 
sublattice A of the mixed-spin Ising model on the decorated Bethe lattice 

m A = o ((^l) _ (^2)) = - ((cTfcl) - (0-fc 2 )) B cthc = ^Bcthc(/3|i?|), (9) 
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which can also be expressed only in terms of the spontaneous magnetization 
WBethe of the corresponding spin- 1/2 Ising model on the undecorated Bethe 
lattice. The other spontaneous sublattice magnetization of the mixed-spin Ising 
model on the decorated Bethe lattice can be descended from the generalized 
Callen-Suzuki identity 35|, |36|, |37[ , which enables to relate it to the previously 
derived spontaneous sublattice magnetization rtiA 



nexp((3Dn 2 ) sinh(/3Jn) 



m B = (S k ) = 2m A 



n=-S 



exp(/3Dn 2 )cosh(/3Jn) 

n=-S 



(10) 



To complete our calculation of spontaneous sublattice magnetizations, one 
still needs the exact closed-form formula for the spontaneous magnetization 
TiBethe of the corresponding spin- 1/2 Ising model on the simple (undecorated) 
Bethe lattice. This quantity can be rigorously found within the framework of 
exact recursion relations [27[. If the undecorated Bethe lattice (see Fig. 1) is 
'cut' at a central site with the spin aki, it will consequently split into q identical 
branches. The partition function of the spin-1/2 Ising model on the simple 
Bethe lattice will take the form 





9 




9n GO. 


+ 


9n vv. 



(11) 



where g n (&ki) is the partition function of the one individual branch and n is 
the total number of generations in the recursively built Bethe lattice. Cutting 
each branch at the subsequent site ak2 allows one to find the recursive relation 
between the partition function of the branch g n (<7ki) with n generations of spins 
and respectively, the partition function of the sub-branch g n _ ifcjyj) with n — 1 
generations of spins 



fl'n(cfci) = y^exp(/3i?g fc io- fc2 ) [g n -l{vk2)] q 



(12) 



The canonical ensemble average of the central spin a^x in the spin-1/2 Ising 
model on the Bethe lattice can readily be calculated from the relation 



^Botho 



1 



Bethe 



ethc 



E 



crfci [g n (cTkx)} q = - 



i n-xi 



1 + xl 



(13) 



which contains the newly defined parameter x n 



g~(-*) 



Even though the 



parameter x n does not have a direct physical meaning, it directly determines 
the overall thermodynamics of the spin-1/2 Ising model on the Bethe lattice in 
the limit n — > oo. With the help of Eq. (TT21) . one may easily derive for x n that 
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enters the magnetization formula (|13[) the following recursion relation 



expf-^Wexpf^W^ 

' ' (14) 



" e X p(M) +ex p(-^)<-l' 

which can be regarded as the equation of state. The equation (fl~4| forms an 
iteration sequence and converges to stable fixed points in the thermodynamic 
limit for non-staggered phases. In our case this situation appears if the effective 
nearest-neighbour coupling of the corresponding spin- 1/2 Ising model on the 
Bethe lattice is ferromagnetic (i.e. R > 0). On the other hand, the recursion 
relation (fT4|) leads to two-cycle double points in the case of staggered antifer- 
romagnetic phase, which emerges in the corresponding spin-1/2 Ising model on 
the Bethe lattice provided that the effective interaction is antiferromagnetic (i.e. 
R < 0) [3^, 4(1 4l|, 42 1 . However, one may take advantage of the bipartite nature 



of the Bethe lattice to overcome this difficulty under the simultaneous change 
of the nearest-neighbour interaction R — > —R and the spins from each second 
generation <n — » (— l)* +1 <7j. Hence, it follows that the staggered magnetization 
|9]) may alternatively be calculated from the relevant exact result (fl~3|) for the 
ferromagnetic spin-1/2 Ising model on the Bethe lattice with the effective in- 
teraction \R\. In this way, our exact calculation of all spontaneous sublattice 
magnetizations of the mixed-spin Ising model on the decorated Bethe lattice 
is completed, because it is now sufficient to substitute Eqs. (TT3l) - (| 1 4 1) into the 
formerly derived expressions ([8|)- (fT0| for the sublattice magnetizations tua, m s A 
and tub- In doing so, one should bear in mind that the effective temperature- 
dependent coupling (3R given by Eqs. ([SJ)— must enter into the recurrence 
relation fT4")) in order to gain the correct exact result for the spontaneous mag- 
netization of the corresponding spin-1/2 Ising model on the simple Bethe lattice 
given by Eq. (fT3|) . 

Finally, let us examine in detail the critical and compensation behaviour of 
the model under investigation. It can be readily understood from the mapping 
relation ([7]) between both partition functions that the mixed-spin Ising model 
on the decorated Bethe lattice may exhibit a critical point only if the corre- 
sponding spin-1/2 Ising model on the undecorated Bethe lattice is at a critical 
point as well. This result is taken to mean that the critical temperature of the 
mixed-spin Ising model on the decorated Bethe lattice can easily be obtained by 
comparing the effective nearest-neighbour coupling of the corresponding spin- 
1/2 Ising model on the simple (undecorated) Bethe lattice with its critical value 
(3 C \R\ = 21n(-3j). The critical condition for the the mixed-spin Ising model on 
the decoratec? Bethe lattice then reads 

s 

exp((3 c Dn 2 ) cosh(^ c Jn) 

e*p(¥)-^ -r^r. W 



exp(&£>n 2 ) 



q-2 



n=-S 
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where (3 C — l/(k-BT c ), T c denotes the critical temperature and the plus (minus) 
sign applies when the effective coupling /3R given by Eqs. ©-([6]) is being 
positive (negative). The sign ambiguity to emerge in the critical condition (|15[) 
comes from the fact that the critical temperature of the spin- 1/2 Ising model 
on the undecorated Bethe lattice remains the same independently of whether 
the effective interaction is being ferromagnetic (R > 0) or antiferromagnetic 
(R < 0). It is worthy to notice, moreover, that the nearest-neighbour spin-1/2 
Ising atoms from the sublattice A should consequently exhibit ferromagnetic 
(antiferromagnetic) spin alignment in the phases delimited by critical phase 
boundaries (|15p with plus (minus) sign. 

Last but not least, the compensation phenomenon may also come into play 
if the total magnetization of the mixed-spin Ising model disappears at tempera- 
ture^) below the critical point. To explore this effect, let us therefore define the 
total magnetization per one spin wit = (2m J 4 + gms)/(l-|-g , 5'), which is normal- 
ized with respect to its saturation value m s = (1 + qS)/{2 + q). In this respect, 
one easily finds the necessary (but not sufficient) condition for an appearance 
of the compensation point (s) 

s 

^ nexp(f3 comp Dn 2 )s'mh(f3 comp Jn) 

1 + 9^ =0, (16) 

22 exp(l3 comp Dn 2 )cosh(l3 com pJn) 

n=-S 

whereas /3 comp — \/(k^T comp ) and T comp labels the compensation temperature. 
It is quite obvious from Eq. (|16|) that the compensation point may possibly arise 
only in the ferrimagnetic model with the antiferromagnetic (J < 0) nearest- 
neighbour interaction between the spin-1/2 and spin-S* atoms. In addition, 
it also turns out from Eq. (1161) that the compensation temperature does not 
apparently depend on a strength of the next-nearest-neighbour interaction J' 
between the spin-1/2 Ising atoms, which merely determines a magnitude of the 
critical temperature according to Eq. (|15j) . 

3. Results and Discussion 

Let us proceed to a discussion of the most interesting results. Even though 
all the results derived in the foregoing section are valid for arbitrary value of the 
coordination number q of the underlying Bethe lattice, regardless of whether the 
interactions J and J' are ferromagnetic or antiferromagnetic and also irrespec- 
tive of the value of the decorating spins, our further analysis will be henceforth 
restricted only to the particular case of the ferrimagnetic mixed spin-1/2 and 
spin-1 Ising model on the decorated Bethe lattice with a quite general coordi- 
nation number q, the antiferromagnetic (J < 0) nearest-neighbour interaction 
and either ferromagnetic (J 7 > 0) or antiferromagnetic (J 7 < 0) next-nearest- 
neighbour interaction. It is noteworthy that the investigated model displays all 



8 



general features of ferrimagnetism found in the mixed spin-(l/2,5) Ising models 
with the integer-valued decorating spins S, while the respective behaviour of the 
mixed spin- (1/2, S 1 ) Ising models with the half-odd-integer decorating spins S is 
much less pronounced as convincingly evidenced in Refs. 0, H, M, E3] and will 
be therefore postponed from the present study. Note furthermore that exact 
results presented hereafter may also be confronted with several known rigor- 
ous results, which have been previously obtained for the mixed spin-1/2 and 
spin-1 Ising ferrimagnets on the decorated planar lattices 0, or Bethe lat- 
tices H, 43, 44 1 . To reduce the number of free parameters, all the interaction 
constants will be hereafter normalized with respect to a strength of the nearest- 
neighbour interaction |J|, which means that the ratio J'/\J\, D/\J\ and k-oT/\ J\ 
will accordingly measure a relative magnitude of the next-nearest-neighbour in- 
teraction, single-ion anisotropy and temperature, respectively. 

Let us begin with the analysis of the ground state. The ground-state phase 
diagram of the mixed spin-1/2 and spin-1 Ising ferrimagnet on the decorated 
Bethe lattice of arbitrary coordination number q is displayed in Fig. [2] In 
agreement with our expectations, the typical ferrimagnetic phase (FI) with 
the antiparallel alignment [cr^Si] = [1/2,-1] between the nearest-neighbour 
spin-1/2 and spin-1 atoms is being the ground state provided that a relative 
strength of the easy-plane (negative) single-ion anisotropy D/\J\ and/or the 
antifcrromagnetic next-nearest-neighbour interaction J'/\J\ is not too strong. 
On assumption that the next-nearest-neighbour interaction is being ferromag- 
netic (J'/\J\ > 0), the decrease in the single-ion anisotropy causes a first-order 
phase transition from FI towards the remarkable ferromagnetic phase (FII) 
with the respective spin alignment [oi,Si] = [1/2,0]. Note that the sponta- 
neous long-range order inherent to FII is merely kept by the ferromagnetic 
next-nearest-neighbour interaction between the spin-1/2 atoms, whereas all the 
spin-1 atoms are non-magnetic as a result of sufficiently strong easy-plane single- 
ion anisotropy (D/\J\ < —1.0) that energetically favours the non-magnetic spin 
state Si — before two magnetic ones Sj = ±1. If one considers the anti- 
ferromagnetic next-nearest-neighbour interaction from the intermediate range 
J'/\J\ £ (—2.0,0.0), the antiferromagnetic phase (AFI) characterized through a 
more complex spin arrangement [<7j, Sf, CTi+i, Si+i] = [1/2, 0; —1/2, 0] arises from 
FI upon lowering the single-ion anisotropy. All the decorating spin-1 atoms re- 
side in AFI the non-magnetic spin state Si — quite similarly as in FII, but 
the spin-1/2 atoms from the nearest-neighbour shells of the decorated Bethe 
lattice align antiparallel with respect to each other on account of the antifer- 
romagnetic nature of the next-nearest-neighbour interaction being responsible 
for the symmetry breaking (period doubling). It should be noted here that 
the emergence of both FII and AFI is restricted to the parameter space with 
the easy-plane character of the single-ion anisotropy D/\J\ < 0. Considering 
the easy-axis single- ion anisotropy D/\J\ > and reinforcing the antiferromag- 
netic next-nearest-neighbour interaction, the mixed-spin system undergoes at 
J'/\J\ = —2.0 a first-order phase transition from FI towards the partially anti- 
ferromagnetically ordered and partially disordered phase (AFII) characterized 
through the spin arrangement [at, Si] <Ji+i, St+i] = [1/2, ±1; —1/2, ±1]. The pe- 
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Figure 2: Ground-state phase diagram of the mixed spin-1/2 and spin-1 Ising model on the 
decorated Bethe lattice in the J'/\J\ — D/\J\ plane. Square brackets [ar,S\ demonstrate the 
relevant spin alignment in each sector of the ground-state phase diagram. For the more 
detailed description of individual phases see the text. 

culiar coexistence of a partial order and a partial disorder in AFII comes from the 
antifcrromagnetic alignment of the next-nearest-neighbouring spin-1/2 atoms, 
which is simultaneously responsible for a unusual spin frustration of the deco- 
rating spin-1 atoms. Hence, it follows that the spin-1 atoms reside in AFII with 
the same probability one of its two magnetic spin states Si = ±1 due to the 
easy-axis character of the single-ion anisotropy D/\J\ > 0. It should be pointed 
out that the established ground-state phase diagram holds for the mixed-spin 
Ising models on the decorated planar lattices as well. It surprisingly appears 
that the characterization of the phase FII has been notably overlooked in several 
previous approximative 0, H| as well as exact 0, [l(| studies on the mixed-spin 
Ising model on decorated planar lattices. In addition, it seems that the partic- 
ular case with the antiferromagnetic next-nearest-neighbour interaction has not 
been dealt with hitherto and hence, theoretical predictions of these two phases 
represent the completely new finding not reported in the literature yet. 

Now, let us focus our attention on a discussion of finite-temperature phase 
diagrams, which are plotted in Figs. [3j6] in the form of critical temperature 
versus single-ion anisotropy dependences for several values of the coordina- 
tion number q and different strengths of the next-nearest-neighbour interac- 
tion J'/\J\. It is noteworthy that all displayed phase boundaries are the lines 
of second-order phase transitions and the ordered (disordered) phases can be 
located below (above) depicted critical boundaries. Fig. [3](a) shows the crit- 
ical temperature as a function of the single-ion anisotropy in an absence of 
the next-nearest-neighbour interaction (i.e. J'/\J\ = 0.0). In this particu- 
lar case, the spontaneously long-range ordered FI represents the ground state 
for D/\J\ > —1.0, while the disordered paramagnetic phase (PP) becomes the 
ground state for D/\ J\ < —1.0. As one can see, the critical line approaches the 
zero-temperature phase boundary between FI and PP with a negative slope for 
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Figure 3: Critical temperature of the mixed spin- 1/2 and spin-1 Ising model on the decorated 
Bethe lattice as a function of the single-ion anisotropy D/\J\ for: (a) J' /\J\ = 0.0 and q=3, 
4, 5 and 6; (b) q = 4 and J'/\J\ < 0. 



the decorated Bethe lattices with the coordination number q < 4, with a posi- 
tive slope for the decorated Bethe lattices with the coordination number q > 4 
and with an infinite gradient for the particular value of the coordination number 
gr = 4. These observations would suggest that reentrant phase transitions can be 
observed in a close vicinity of the ground-state boundary between FI and PP for 
the decorated Bethe lattices with a sufficiently high coordination number q > 4, 
whereas the reentrance becomes the more pronounced, the higher the coordi- 
nation number q of the underlying Bethe lattice is. The displayed exact results 
may be contrasted with the approximative results reported by Kaneyoshi 0, Q , 
which imply either more 0] or less [§| robust reentrant region for the decorated 
square lattice. On the other hand, the present results are in accordance with 
the previously reported exact results for mixed-spin Ising model on the deco- 
rated planar lattices, among which the decorated square [iiEll, honeycomb and 
kagome lattice [H} do not show reentrance in contrast to the decorated trian- 
gular lattice with a pronounced reentrant region 15;]. Note furthermore that 
the investigated model cannot exhibit the compensation phenomenon when the 
next-nearest-neighbour interaction is absent. 

Next, we will pay our attention to the effect of the antiferromagnetic next- 
nearest-neighbour interaction on the critical behaviour of the mixed-spin Ising 
model on the decorated Bethe lattice with the coordination number q — 4. It can 
be clearly seen from Fig. Elb) that the critical temperature of FI monotonon- 
ically decreases with decreasing the single-ion anisotropy until it reaches the 
zero-temperature phase boundary between FI and AFI (see the right wing of 
each displayed critical lines in Fig. EJb)). Afterwards, the new critical line 
develops for a more negative values of the single-ion anisotropy, whereas the op- 
posite trend is observed; the critical temperature monotonically increases with 
decreasing the single-ion anisotropy until it tends towards its maximum value in 
the limit D/\ J\ — > —oo. The left wing of each critical line obviously corresponds 
to AFI phase with the zero spontaneous magnetization but the non-zero stag- 
gered magnetization of the spin-1/2 Ising atoms (see the subsequent discussion 
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(a) DI\J\ (b) DI\.J\ 



Figure 4: Critical temperature of the mixed spin-1/2 and spin-1 Ising model on the decorated 
Bethe lattice with the coordination number q = 4 as a function of the single-ion anisotropy 
D/\ J\ for several values of the ferromagnetic next-nearest-neighbour interaction J'/\J\ > 0. 
Broken curve shows the line of compensations points and Fig. HJb) is a detailed plot from the 
region near the ground-state phase boundary between FI and FII. 

for more details). It should be noticed that the antiferromagnetic next-nearest- 
neighbour interaction can neither induce the reentrant phase transitions, nor 
the compensation phenomenon for the decorated Bethe lattices with the coor- 
dination number q < 5. 

For completeness, the finite-temperature phase diagram of the mixed spin- 
1/2 and spin-1 Ising model on the decorated Bethe lattice with the coordination 
number q = 4 is plotted in Fig. [4] for different strengths of the ferromagnetic 
next-nearest-neighbour interaction. The critical temperature generally exhibits 
a smooth monotonous decline with decreasing the single-ion anisotropy and at 
first sight, there is no clear evidence of the phase transition between FI and FII 
phases in the displayed critical lines. It will be later demonstrated, however, 
that the sublattice magnetization of the spin-1 atoms approaches in the zero- 
temperature limit its saturation value only if D/\ J\ > —1.0, while it starts from 
zero in the reverse case. With this in mind, one may conclude that the critical 
line of FI continuously merges at some finite temperature with the critical line of 
FII at the specific value of the single- ion anisotropy D/\ J\ = —1.0, whereas the 
former FI phase is located below the right part (D/\J\ > —1.0) and the latter 
FII phase below the left part (D/\J\ < —1.0) of the overall critical line. Last but 
not least, it is worth mentioning that a presence of the weak (strong) ferromag- 
netic next-nearest-neighbour interaction may lead to an existence of one (two) 
compensation points in the temperature dependence of the total magnetiza- 
tion as it can be clearly seen from the respective behaviour of the broken curve, 
which displays the numerical solution of the necessary condition (| 16[) for appear- 
ance of the compensation points. It should be remembered that the necessary 
condition for appearance of compensation points does not depend according to 
Eq. (fT6|l on a relative strength of the next-nearest-neighbour interaction J'/\J\. 
From this point of view, the line of compensation temperatures starts at the 
ground-state boundary between FI and FII phases at D/\J\ = —1.0 and it 
ends up at a respective crossing point between the necessary condition (|16|) and 
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Figure 5: Critical temperature of the mixed spin- 1/2 and spin-1 Ising model on the decorated 
Bethe lattice with the coordination number q = 6 as a function of the single-ion anisotropy 
D/\J\ for several values of the antiferromagnetic next-nearest- neighbour interaction J'/\J\ < 
0. Broken curve shows the line of compensations points and Fig. [5jb) is a detailed plot from 
the observed reentrant region near the ground-state phase boundary between FI and AFI. 

the relevant critical line depending basically on a relative strength of the next- 
nearest-neighbour interaction (above this line the spin system is disordered and 
there is no spontaneous ordering). On the other hand, it also turns out that 
the ferromagnetic next-nearest-neighbour interaction cannot induce reentrant 
phase transitions for the decorated Bethe lattices with the coordination number 
q < 5. 

To provide a deeper insight into how the coordination number affects the 
overall critical behaviour, the critical temperature of the mixed spin-1/2 and 
spin-1 Ising model on the decorated Bethe lattice with the higher coordination 
number q = 6 is plotted against the single-ion anisotropy in Figs. [5JS1 First, let 
us briefly comment on the most significant changes of the relevant critical lines 
when considering the influence of the antiferromagnetic next-nearest-neighbour 
interaction (Fig. [5J. In the absence of the next-nearest-neighbour interaction, 
the investigated spin system exhibits a relatively wide reentrant region for the 
single-ion anisotropy close to but slightly below the ground-state boundary be- 
tween FI and PP at D/\J\ = —1.0. In this special case, the spin system is 
disordered in the ground state, then it orders at a lower critical temperature 
to display the spin order inherent to FI, which finally disappears at the upper 
critical temperature (PP-FI-PP reentrance). On the other hand, there are quite 
fundamental differences in the respective reentrant behaviour of the spin system 
with the antiferromagnetic next-nearest-neighbour interaction. First of all, the 
spin system exhibits a more striking reentrance with three consecutive critical 
points when starting from the spontaneously long-range ordered ground state, 
which either corresponds to AFI or FI depending on an interplay between the 
single-ion anisotropy and the next-nearest-neighbour interaction. The spon- 
taneous ordering of AFI (or FI) disappears at the first critical temperature, 
then there appears PP in a relatively narrow region in between the first and 
second critical temperature, at which the spontaneous order of FI occurs (or 



13 



-3 -2 -1 1 -1.50 -1.25 -1.00 -0.75 

(a) DI\J\ (b) DI\J\ 

Figure 6: Critical temperature of the mixed spin- 1/2 and spin-1 Ising model on the decorated 
Bethe lattice with the coordination number q = 4 as a function of the single-ion anisotropy 
D/\ J\ for several values of the ferromagnetic next-nearest-neighbour interaction J' /\J\ > 0. 
Black broken curve shows the line of compensations points, whereas Fig. [4fb) depicts a detail 
from the region, where reentrant phase transitions can be observed. 

re-appears) and finally vanishes at the third critical temperature. The observed 
reentrant phase transitions can be thus characterized by two following sequences 
of second-order phase transformations AFI-PP-FI-PP and FI-PP-FI-PP. It can 
be also clearly seen from Fig. [5jb) that the reentrant region moves towards a 
more positive values of the single- ion anisotropy (together with the ground-state 
boundary between FI and AFI) and it gradually vanishes upon strengthening 
the antiferromagnetic next-nearest-neighbour interaction. 

Finally, the last couple of phase diagrams from Fig. [6] illustrates the effect of 
the ferromagnetic next-nearest-neighbour interaction on the critical behaviour 
of the mixed spin- 1/2 and spin-1 Ising model on the decorated Bethe lattice 
with the coordination number q — 6. As one can see, the character of reen- 
trant phase transitions substantially changes by taking the ferromagnetic next- 
nearest-neighbour interaction into consideration even although one still detects 
the reentrance with three consecutive critical points. However, the observed 
reentrant phase transition is basically different in that both low-temperature as 
well as high-temperature spontaneously long-range ordered phase is always FII. 
According to this, the spin system undergoes the following sequence of second- 
order phase transitions FII-PP-FII-PP when starting from the lowest and end- 
ing up at the highest temperature. To summarize, the only common feature of 
all three kinds of reentrant phase transitions FII-PP-FII-PP, FI-PP-FI-PP and 
AFI-PP-FI-PP is that the reentrance is quite sensitive and easily diminishes by 
changing a relative strength of the next-nearest-neighbour interaction. 

To provide an independent check of the aforedescribed critical behaviour, 
let us proceed to a discussion of thermal dependences of the order parame- 
ters as depicted in Figs. [TJIHl First, let us confirm the phase transition from 
FI towards FII by considering the constant value of the ferromagnetic next- 
nearest-neighbour interaction J' /\J\ = 3.0 and changing a relative strength of 
the single-ion anisotropy. Fig. [Zta)-(b) illustrate some typical temperature 
dependences of the total and sublattice magnetizations when FI is being the 
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(c) k B T/\J\ (d) k B T/\J\ 



Figure 7: Thermal variations of the total and sublattice magnetizations for the mixed spin- 
1/2 and spin-1 Ising model on the decorated Bethe lattice with the coordination number 
q = 4, the fixed value of the ferromagnetic next-nearest-neighbour interaction J'/\J\ = 3.0 
and four different values of the single-ion anisotropy: (a) D/\ J\ = 1.0; (b) D/\ J\ = —0.9; (c) 
D/\J\ = -1.0; (b) D/\J\ = -1.5. 



ground state. The resultant magnetization exhibits for D/\J\ = 1.0 and —0.9 
the standard R-type dependence and the slightly deformed N-type dependence 
with one compensation point, respectively, whose unusual shape comes from 
a steeper temperature-induced decrease in the sublattice magnetization ms of 
the decorating spin-1 atoms. Furthermore, Fig. EJc) shows temperature varia- 
tions of the total and sublattice magnetizations exactly at the phase boundary 
between FI and FII. In agreement with this statement, the sublattice magne- 
tization rriB tends towards the intermediate value 0.5 when approaching zero 
temperature, which indicates a coexistence of FI and FII in the ground state. Fi- 
nally, thermal dependences of the total and sublattice magnetizations displayed 
in Fig. [7jd) provide an exact evidence of FII. As a matter of fact, the sublat- 
tice magnetization tub vanishes in the zero-temperature limit and its anoma- 
lous temperature-induced increase can be attributed to spin excitations from 
the non-magnetic spin state Si — preferably towards the excited spin state 
Si = — 1. In the consequence of that, the resultant magnetization might display 
the striking W-type dependence with two successive compensation points not 
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(a) k B T/\J\ (b) k B T/\J\ 



Figure 8: Thermal variations of the total, sublattice and staggered magnetizations for the 
mixed spin- 1/2 and spin-1 Ising model on the decorated Bethe lattice with the coordination 
number q = 4, the fixed value of the antiferromagnetic next-nearest-neighbour interaction 
J'/| J| = —0.2 and two different values of the single-ion anisotropy: (a) D/\J\ = —0.8; (b) 
D/\J\ = -1.0. 



included in the standard Neel theory of ferrimagnetism (for the characterization 



of temperature dependences of the total magnetization see Refs. [4a, |46(). The 
former compensation point appears owing to a higher number of the decorating 
spins (the overall contribution of the sublattice magnetization raj to the to- 
tal magnetization consequently exceeds the one of the sublattice magnetization 
m,A provided that tub becomes sufficiently high), while the latter compensation 
point emerges as a result of more subtle character of the decorating spins with 
respect to thermal fluctuations (the sublattice magnetization mg tends steeper 
to zero than the sublattice magnetization rriA in a vicinity of the critical point). 

Next, let us turn to thermal variations of the total and sublattice magnetiza- 
tions of the mixed spin-1/2 and spin-1 Ising model on the decorated Bethe lattice 
with the coordination number q — 4 and the antiferromagnetic next-nearest- 
neighbour interaction J'/\J\ = —0.2. It has been demonstrated that the an- 
tiferromagnetic next-nearest-neighbour interaction competes with the nearest- 
neighbour interaction and hence, the weaker easy-plane single-ion anisotropy is 
thus generally needed in order to destroy the spontaneous order of FI. However, 
one still detects temperature dependences of the total and sublattice magneti- 
zations quite typical for FI whenever the single-ion anisotropy is selected above 
its critical value Df,/\J\ = —1 — J'/ (2| J|) (see Fig. EJa)). Contrary to this, both 
spontaneous sublattice magnetizations equal zero for any single-ion anisotropy 
below this boundary value D < Dt and the system surprisingly exhibits the 
unusual antiferromagnetic long-range order AFI with the non-zero staggered 
magnetization of the spin-1/2 atoms as the relevant order parameter (see Fig. 
E^b)). Hence, it actually turns out that the spin-1/2 Ising atoms from nearest- 
neighbour shells of the decorated Bethe lattice align antiparallel with respect to 
each other due to the antiferromagnetic next-nearest-neighbour interaction and 
owing to this fact, the decorating spin-1 atoms become frustrated and all tend 
towards their non-magnetic spin state Si = 0. 
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(a) k B 77 \J\ (b) KTI\J\ 



Figure 9: Thermal variations of the total, sublattice and staggered magnetizations for the 
mixed spin- 1/2 and spin-1 Ising model on the decorated Bethe lattice with the coordination 
number q = 6, the fixed value of the antiferromagnetic next-nearest-neighbour interaction 
J'/|J| = —0.2 and two different values of the single-ion anisotropy: (a) D/\J\ = —0.89; (b) 
D/\J\ = -0.91. 



For completeness, let us also confirm an existence of the reentrant phase 
transitions by analyzing thermal dependences of the spontaneous magnetization 
of the mixed spin-1/2 and spin-1 Ising model on the decorated Bethe lattice with 
the higher coordination number q = 6. For this purpose, we depict in Fig. [9] 
temperature variations of the total and sublattice magnetizations for one specific 
value of the antiferromagnetic next-nearest-neighbour interaction J'/\J\ = —0.2 
and two different values of the single-ion anisotropy. In the former case with 
D/\J\ — —0.89, one observes three successive reentrant phase transitions just 
in between two phases FI and PP in the following order FI-PP-FI-PP. On the 
other hand, another type of reentrance can be detected in the latter case with 
D/\J\ — —0.89, which reveals three consecutive phase transitions between AFI, 
FI and PP in the following order AFI-PP-FI-PP. It should be stressed that 
the latter reentrance would be misinterpreted as the reentrant phase transition 
with only two consecutive critical points if the staggered magnetization would 
not be taken into account as the possible order parameter along with the uniform 
magnetization. 



4. Conclusion 

In this work, we have exactly calculated magnetic properties of the mixed 
spin-1/2 and spin-S" Ising model on the decorated Bethe lattice with a quite 
general coordination number q when combining the decoration-iteration trans- 
formation with the exact recursion relations. The model Hamiltonian includes 
the nearest-neighbour interaction between the spin-1/2 and spin-S" atoms, the 
next- nearest- neighbour interaction between the spin-1/2 atoms and the uniaxial 
single-ion anisotropy. By making use of the decoration-iteration transformation, 
we have established a precise mapping equivalence between the mixed spin-1/2 
and spin-S 1 Ising model on the decorated Bethe lattice and its corresponding 
spin-1/2 Ising model on the simple (undecorated) Bethe lattice subsequently 
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treated with the aid of exact recursion relations. Furthermore, the rigorous 
mapping theorems and exact spin identities have been employed in order to get 
the precise results for both spontaneous sublattice magnetizations as well as 
staggered magnetization. 

Our main attention was focused on magnetic properties of the mixed spin- 
1/2 and spin-1 Ising ferrimagnet, which exhibits remarkably diverse critical 
behaviour and compensation phenomena. The four different phases have been 
found in the ground-state phase diagram, two of them have the uniform mag- 
netization and another two staggered magnetization as the order parameter. 
In addition, it has been demonstrated that the investigated mixed-spin sys- 
tem may exhibit either reentrant phase transitions with two consecutive critical 
points (for zero next-nearest-neighbour interaction) or three different types of 
reentrant phase transitions with three successive critical points (for non-zero 
next-nearest-neighbour interaction). It is worthy to recall that the reentrance 
may appear just for a sufficiently high coordination number q > 4 of the under- 
lying Bethe lattice and this observation seems to be of a quite general validity as 
it has been already proved in the Ising-Heisenberg models on the diamond-like 



decorated Bethe lattices [47|, |48[ . Finally, it has been also illustrated that the 
next-nearest-neighbour interaction may cause a presence of one or two compen- 
sation points in the temperature dependence of the total magnetization. 

Let us conclude our study by mentioning that the present formalism can 
be rather straightforwardly extended for the mixed spin- 1/2 and spin- 5 Ising 
model on the decorated Bethe lattices in a presence of the external magnetic field 
and this exactly soluble model might display a quite intriguing magnetization 
process including several non-trivial magnetization plateaus. In this direction 
will continue our further efforts. 
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